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ABSTRACT 

The current status of theoretical QCD calculations and experimental measurements of the 
Gottfried sum rule are discussed. The interesting from our point of view opened problems 
are summarised. Among them is the task of estimating the measure of light-quark flavour 
asymmetry in possible future experiments. 
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1 Introduction. 



Studies of the Gottfried sum rule of charged lepton-nucleon deep- inelastic scattering [T], 
namely 

I G (Q 2 , 0, 1) = f 1 - \f?(x, Q 2 ) - Ff(x, Q 2 )] (1) 

JO XL 

can provide the important information on the possible existence of a light antiquark flavour 
asymmetry in the nucleon sea. Indeed, the NMC collaboration determination [2] demon- 
strated that its experimental value 

l£ MC (4 GeV 2 , 0,1) = 0.235 ±0.026 , (2) 

is significantly lower than the quark-parton flavour-symmetric prediction 

/g(Q 2 ,0,1) = 1 . (3 ) 

This deviation is associated with the existence of a non-zero integrated light-quark flavour 
asymmetry defined as 

FA(Q 2 , 0, 1) = [\d(x, Q 2 ) - u(x, Q 2 )]dx . (4) 
Jo 

In spite of the existence of detailed reviews on the subject |HJ E] we think that it is worth 
while to return to the consideration of the current status of knowledge of different aspects 
related to this sum rule. 

In this report, based in part on the recent work of Ref. 0, the contributions of QCD 
corrections and higher-twist effects to the Gottfried sum rule are discussed first in the 
case of a flavour-symmetric sea. Next, the results of its experimental determination are 
summarised. Then we briefly outline various possibilities for determinating the integral 
FA(Q 2 , 0, 1) from previos, present and future data. 



2 QCD predictions. 



Let us start by defining the arbitrary non-singlet (NS) Mellin moment of the difference of 
F 2 structure functions of charged lepton-proton and charged lepton-nucleon deep-inelastic 
scattering (DIS): 



n-2 



x 



F^(x,Q 2 )-Ff(x,Q 



dx 



(5) 



The moment with n = 1, namely the Gottfried sum rule, can be expressed as 



/ G (Q 2 ,o,i) 
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(u v (x, Q 2 ) - d v {x, Q 2 )) + -{u(x, Q 2 ) - d{x, Q 2 )) 



dx 



---FA(g 2 ,o,i) 



(6) 



where u v (x, Q 2 ) = u(x, Q 2 ) —u(x, Q 2 ) and d v (x, Q 2 ) = d(x, Q 2 ) — d(x, Q 2 ) are the valence- 
quark distributions and the measure of flavour asymmetry is related to the difference of 
the sea-quark distributions u(x, Q 2 ) and d(x, Q 2 ) via Eq.(JH). 



2.1 Perturbative contributions. 



Consider first the case when the sea is flavour- symmetric. In zeroth order of perturba- 
tion theory the quark-parton result of Eq.(j3J) is reproduced. However, the quark-gluon 
interactions generate non-zero corrections to Iq, defined as 



J G (Q 2 , 0,1) = AD(a s )C(a s 



(7) 



The anomalous dimension term is related to the anomalous dimension function of the first 



moment 7 n (a s ) and to the QCD /3-function, namely : 
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The given expansion in a s can be obtained after taking into account that 7q =1 = and 
setting 5 = 0. 

The calculations of 7™ =1 [3 EI give the following result 

-4(C 2 - C F C A )[13 + 8C(3) - 2tt 2 



7P 1 



+2.557 



The coefficients of the QCD /3-function are well-known: 
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(9) 

(10) 
(11) 



Here and below Cp = 4/3 and Ca = 3 and / is the number of active flavours. 
The general order a 2 - expression for C(a s ) can be written down as 



C(a s 



1 + C\ 



n=l 



+ a 



n=l 



(12) 



where C™ =1 = |8 3 . The coefficient C% =1 was evaluated only recently [Sj by means of 
numerical integration of the complicated x-dependence of the two-loop contributions to 
the coefficient functions of DGLAP equation for DIS structure functions, calculated in 
Ref. p. The expression obtained in Ref. [S] is: 



cr 1 



(3.695C| - IM7C F C A ) = -0.821 
Collecting now all known QCD corrections to Eq.(7) we found [H|: 



j g (q 2 ,o,i) /=3 
/ G (g 2 ,o,i) /=4 
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(13) 

(14) 
(15) 



where a s = a s (Q 2 ) and the three-loop anomalous dimension term 72 =1 is still unknown. In 
Ref. [5] it was estimated using the feature observed in Ref. [12] that the n-dependence of the 
ratio 7^/72 , obtained from three- loop terms of the the anomalous dimension functions of 
even moments for charged lepton- nucleon DIS, calculated in Ref.|l(Jj. and of odd moments 
of vN DIS, calculated in Ref. ^T], can be fixed by similar approximate relation. Taking 
these estimates into account we got [3] 



/ G (Q 2 ,0,l) /= 



1 + 0.0355 



0.811(^ 

V 7T 



2i 



(16) 



/ G (Q 2 ,o,i) /= 
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(17) 



where the a 2 contribution is dominated by the numerical value of the coefficient C% =1 
from Eq. (|13j) . Thus we convinced ourselves that the perturbative QCD corrections to 
the Gottfried sum rule are really small and cannot be responsible for violation of the 
flavour-symmetric prediction from the experimental value of Eq.(2). 



2.2 Higher- twist terms. 

The possibility that the higher-twist effects in the Gottfried sum rule might be sizeable 
was discussed in Ref. [13]. It was argued that the next-to- leading sets of parton distribu- 
tions, namely GRV94 [H], MRST98 [J5] and CTEQ5 [To], failed to describe the existing 
experimental F^-F^ data below Q 2 < 7 GeV 2 [13]. From the point of view of the authors 
of Ref. ^3] this might be associated with substantial higher-twist corrections, which in 
part are responsible for the deviation of the Gottfried sum rule result from its NMC value. 
However, definite results of the fits to Ff — F% data, performed in Ref. ^7] with the help 
of the most recent Alekhin PDF set of Ref. ^B] (A02), indicate that that the conclusions 
of Ref. [T3] are too optimistic. Indeed, the authors of Ref. [T7J demonstrated that for 
the second NS moment of Eq. (5) the numerator of the l/Q 2 twist-4 correction is rather 
small, namely E n F ~ v = -0.0058 ± 0.0069 GeV 2 . In view of this we expect that the twist-4 
contribution to the first moment, namely to J G , will be small also. Moreover, in spite of 
the fact that the fits of Ref. ^7] reveal a definite discrepancy between x-dependence found 
for Hp~ p (x) and the predictions of the infrared renormalon (IRR) model (for a review see 
Ref. fH]), we think that the latter method still might give order-of-magnitude estimates of 
the higher-twist contributions. Note, that the IRR model is based in part on summations 
of the chain of fermion loop insertions to the gluon propagator and is thus related to the 
large /-expansion of the coefficient functions. For the polarised Bjorken sum rule these 
studies were made in Ref. |20j (for more recent discussions see Ref. [21j). In the case of 
Ig the flavour dependence does not manifest itself up to a^-corrections (see Eq. (12) ). 
Thus, we conclude, that in comparison with polarised Bjorken sum rule, the IRR model 
corrections (and therefore higher-twist effects) to the Gottfried sum rule will be damped 
by the additional factor a s /7i. In view of this we think that the Gottfried sum rule cannot 
receive substantial higher twist contributions, though the explicit demonstration of the 
validity of this statement is still missing. 



3 Experimental situation. 



The experimental determinations of the Gottfried sum rule have a rather long history, 
summarised in the reviews of Refs. [3J E] . In fact what is really evaluated from the experi- 
mental data is the integral 



I(Q 2 ,x r , 







m.in 



F^( Xl Q 2 )-F^(x,Q' 



In a more detailed analysis F^^x, Q 2 ) — F!f l {x, Q 2 ) should be extrapolated to the unmea- 
sured regions and since F^ix, Q 2 ) is extracted from DIS on nuclei targets, nuclear effects 
should be also taken into account. However, all four experimental groups working on the 
direct determiantion of the Gottfried sum rule from their experimental data, were not able 
to achieve ideal results. Indeed, the main source of experimental uncertainty results from 
the extrapolations of the experimental data from x m i n to 0. Moreover, the mean Q 2 in 
the data are dependent on x and it is sometimes difficult to fix typical Q 2 value of the 
Gottfried sum rule (see Table 1.) 



Group 


Q 2 (GeV 2 ) 








SLAC 22 


0.1-20 


0.02 


0.82 


0.200±0.040 


EMC [23^ 


10-90 


0.02 


0.8 


0.197±0.011(st.)±0.083 (sys.) 


BCDMS [24 


20 


0.06 


0.8 


0.197±0.011 (st.)±0.036(sys.) 


NMC [2 j 


4 


0.004 


0.8 


0.221±0.008 (st.)±0.019(sys.) 



Table 1. The existing experimental data for the integral of Eq. (|18|). 

In spite of the fact that already the results of Ref. [22] inspired discussions of the 
possibility that the theoretical expression for Eq.(j3J) might be violated, its further deter- 
minations by EMC collaboration j2H], namely 

J™ C (Q 2 =?, 0,1) = 0.23518;$ (19) 

within experimental error bars did not in fact demonstrate an obvious deviation from the 
quark-parton model prediction 1/3 (note, that in Eq.(19) Q 2 -value was not determined). 
A similar conclusion also applies to the analysis of the BCDMS data in Ref. j21]. Indeed, it 
suffers from large uncertainties at x < 0.06, which vary from 0.07 to 0.22. Thus the results 
of NMC collaboration of Eq. (J2J) turned out to be extremely important for understanding 
that flavour-asymmetry of antiquark distributions in the nucleon really exist in nature. 
The precision of their data even allow one to extract the value of integrated light-quark 
flavour asymmetry, defined by Eq.(4) [2]: 

FA NMC (4 GeV 2 , 0,1) = 0.147 ±0.039 . (20) 

However, even the members of NMC collaboration were not able to take into account all 
effects, typical for DIS. Indeed, the nuclear corrections, such as the Fermi motion, were 
neglected by them. In view of this, it became rather important to get an independent 



experimental extraction of FA(Q 2 , 0, 1). Quite recently this was done in an analysis of 
the data for Drell-Yan production in proton-proton and proton-deuteron scattering by the 
members of E866 collaboration. Integrating d(x, Q 2 ) —u(x, Q 2 ) over the measured x-region 
they obtained |25] : 

FA E866 (54 GeV 2 , 0.015, 0.35) = 0.0803 ±0.011 . (21) 

Extrapolation this integral to the unmeasured region < x < 0.015 and assuming that 
the contribution for x > 0.35 is negligible, the members of E866 collaboration found that 

FA E866 (54 GeV 2 , 0, 1) = 0.118 ± 0.015 (22) 

(see Ref. |25j). Within existing error-bars this result turned out to be in agreement with 
Eq. (20), which is the NMC value of this integral, extracted at 4 GeV 2 . In view of this it 
is possible to conclude that the value of the integrated light-quark flavour asymmetry is 
almost independent of Q 2 over a wide range of the momentum transfer. This demonstrates 
in part its non-perturbative origin. 



4 Possible future prospects. 

In order to calculate the characteristics of a light-quark flavour asymmetry a number of 
non-perturbative models have been successfully used. Among them are the meson-cloud 
model, instanton model, chiral-quark soliton models and others (for reviews see Refs. 
[3J H]). However, to make more definite conclusions on the predictive power of these 
approaches more detailed experimental knowledge about the behaviour of the rate d/u at 
different values of x is needed. 

The new 120 GeV Fermilab Main Injector should allow one to extend Drell-Yan mea- 
surements of d/u to the region of 0.02 < x < 0.7. Moreover, the studies of the CEBAF 
data for F® al large x can give the chance to perform more detailed combined fits of all 
available DIS data. The extraction of a light-quark flavour asymmetry from the fits to 
the data for F® require the detailed treatment of nuclear effects, say in the manner of 
the work of Ref. [12] . Clearly, these measurements might be important for more detailed 
determinations of the effects of flavour-asymmetry in various sets of parton distribution 
functions, which at present differ in the CITEQ6M, MRST2001C and A02 sets (for their 
comparison at Q 2 = (100 GeV) 2 see Ref. j2H])- As a test of their current predictive power 
it can be rather useful to use them for calculations of the integral FA(Q 2 ,0, 1), as was 
done in Ref. [I]. 

Another possibility is to study light-quark flavour asymmetry in the process of uN 
DIS at the possible future Neutrino Factories. Indeed, as was noticed in Ref. [S] in the 
quark-parton model the d — u difference can be related to the vN DIS SFs as 



d{x) — u(x) = — 



F?{x) - F? p (x) 



r p (x) - F7(x) 



(23) 



where the s and c distributions are neglected. At the NLO of perturbative QCD the analog 
of Eq.((2SI) was also derived [Hj. Since it is known, that at the Neutrino Factories it will 



be possible to extract from the cross-sections of vN DIS F\ and F3 structure functions 
separately |2Zj, the more precise uN DIS data might be useful for additional estimates of 
the size of the light-quark flavour asymmetry. 
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